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I. INTRODUCTION 

The Standard Model of particle physics (SM) contains 
a large number of free parameters which can only be de- 
termined by experiment. Precision measurements of all 
of these parameters are essential for probing the validity 
range of the model by comparing many other precision 
measurements to SM calculations. Three of the param- 
eters, the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
element \V c b\ [1, 2] and the heavy quark masses nib and 
m c , can be related via Operator Product Expansions 
(OPE) to moments and rates of inclusive distributions 
in semileptonic B meson decays, B — > X c t~T> [3], and 
rare B-meson decays, B — > X s j, where X c and X s de- 
note the hadronic systems with charm and strangeness in 
these final states, respectively. The quantities \V c b\, m-b, 
m Cl and nonperturbative parameters describing effects of 
the strong interaction can be determined from the mea- 
sured rates and moments using expansions in 1/mj and 
the strong coupling constant a s with reliable uncertainty 
estimates. 

Various measurements of moments of the hadronic- 
mass [4-8] and lepton-energy [7, 9,10] spectra in inclusive 
semileptonic decays B — > X c l~v have already been used 
for determinations of \V c b\, nib, m c , and of four strong- 
interaction parameters /Lt^du), fiQ^/j,), p D (p), and p\g{p)- 
The parameters p^(p) and f%(fJ>) are the expectation 
values of the kinetic and chromomagnetic dimension-five 
operators, respectively, and appear at 0(l/m 2 ) in the 
expansion. The parameters p 3 D (p) and p 3 LS (p) are the ex- 
pectation values of the Darwin and spin-orbit dimension- 
six operators, respectively, and appear at 0(l/m 3 ) in the 
expansion [11]. Here, p denotes the Wilson factorization 
scale that separates effects from long- and short-distance 
dynamics. 

Combined fits to the B — > XJrv moments and mo- 
ments of the photon-energy spectrum in B — > X s ^j de- 
cays [12-16] in the context of Heavy Quark Expansions 
(HQE) lead to |V c6 | = (41.96 ± 0.23 ± 0.69) x lO" 3 and 
m b = (4.590 ± 0.025 ± 0.030) GeV/c 2 in the kinetic-mass 
scheme [17] and \V cb \ = (41.78 ± 0.30 ± 0.08) x 10~ 3 and 
m b = (4.701 ± 0.030) GeV/c 2 in the IS scheme [18]. The 
Belle Collaboration has presented similar results in [19]. 

While lepton-energy moments are known with good ac- 
curacy, the precision of the hadronic-mass and photon- 
energy moments is limited by statistics. Therefore, we 
present a new measurement of the hadronic-mass mo- 
ments (m x ) with k = 1, . . . , 6 based on a larger dataset 
than previously used [5]. We also present the first mea- 
surement of the combined hadronic mass-and-energy mo- 
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ments (n x ) with k — 2, 4, 6 as proposed by Gambino and 
Uraltsev [20]. The combined moments (n x ) use the mass 
mx and the energy Ex of the X c system in the B meson 
rest frame of B — ► X C £~V decays, 

n\ = m 2 x c 4 - 2AE X + i 2 , (1) 

with a constant A, here fixed to be 0.65 GeV as pro- 
posed in [20] . They are expected to allow a more reliable 
extraction of the higher-order nonperturbative HQE pa- 
rameters and thus to increase the precision on the ex- 
traction of \V c b\ and the quark masses rrib and m c . All 
moments are determined for different values of the mini- 
mum energy of the charged lepton. 

We update our previous measurement of lepton-energy 
moments [9] using branching fraction measurements for 
background decays in [21] and improving the evaluation 
of systematic uncertainties. 

Finally, we perform a combined fit to the hadronic- 
mass moments, moments of the lepton-energy spectrum, 
and moments of the photon-energy spectrum in decays 
B — > X s j. The fit determines \V c b\, the quark masses 
nib an< 4 m ci the total semileptonic branching fraction 
B{B — > X c l~v\ and the dominant nonperturbative HQE 
parameters p 2 , pfj, pp, and p\ s . An alternative fit to 
the moments of n x , of the lepton-energy, and of the pho- 
ton energy in B — > X s "f, leads to essentially the same 
results. 



II. BABAR DETECTOR AND DATASETS 

The work is based on data collected with the BABAR 
experiment [22] at the PEP-II asymmetric-energy e + e~ 
storage rings [23] at the SLAC National Accelerator Lab- 
oratory. 

The BABAR tracking system used for charged particle 
and vertex reconstruction has two main components: a 
silicon vertex tracker (SVT) and a drift chamber (DCH) , 
both operating within a 1.5-T magnetic field of a super- 
conducting solenoid. The transverse momentum resolu- 
tion is 0.47 % at 1 GeV/c. Photons are identified in an 
electromagnetic calorimeter (EMC) surrounding a detec- 
tor of internally reflected Chcrcnkov light (DIRC), which 
associates Chcrcnkov photons with tracks for particle 
identification (PID). The energy of photons is measured 
with a resolution of 3 % at 1 GeV. Muon candidates are 
identified with the use of the instrumented flux return 
(IFR) of the solenoid. The tracking system, EMC, and 
IFR cover the full azimuthal range and the polar-angle 
range 0.3 < 8 < 2.7 rad in the laboratory frame, cor- 
responding to a coverage of approximately 90% in the 
center-of-mass (cm.) frame, where 9 is the polar angle 
with respect to the electron direction. The DIRC fidu- 
cial volume corresponds to a cm. frame coverage of about 
84%. 

The data sample for the hadronic moments mea- 
surements consists of about 210 fb _1 , corresponding to 
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232 x 10 6 decays T(AS) — > BB. Our previous measure- 
ment of the lepton-energy moments, which is updated 
in this paper, was based on a data sample of about 
51 x 10 6 Y(AS) — * BB decays. This corresponds to an 
integrated luminosity of 47 fb _1 on the T(4S) resonance. 
In addition, about 9 fb _1 of data recorded at an energy 
40 MeV below the resonance (off-resonance) was used in 
the lepton-energy moments measurement for the subtrac- 
tion of background not originating from the T(45). 

We use Monte Carlo (MC) simulated events to deter- 
mine background distributions and to correct for detector 
acceptance and resolution effects. Simulated B-meson 
decays are generated using EvtGen [24]. The simulation 
of the BABAR detector is realized with GEANT4 [25] and 
final state radiation (FSR) is modeled using the PHOTOS 
code [26]. 

In the simulation of semileptonic decays B — > X C £~V 
we use the branching fractions listed in Table I. For 
the dominant decay B — > D*£~v we use a parameter- 
ization of form factors, based on heavy quark effective 
theory (HQET) [27-29]. Its differential rate is described 
by three helicity amplitudes which are expressed by the 
three parameters p 2 , i?i , and i?2 • We choose the follow- 
ing values measured in [30]: R x = 1.18 ± 0.30 ± 0.12, 
R 2 = 0.71 ± 0.22 ± 0.07, and p 2 = 0.91 ± 0.15 ± 0.06. 
The quoted errors reflect statistical and systematic un- 
certainties. For decays B — ► Di~v and for decays to 
the higher-mass states D\, D[, Dq, and DZ we use the 
ISGW2 model [31]. For the decays B -► D^nlV, we use 
the prescription by Goity and Roberts [32] . 

III. RECONSTRUCTION OF SEMILEPTONIC 
DECAYS FOR THE MEASUREMENT OF 
HADRONIC MOMENTS 

The event selection and reconstruction for the 
hadronic-mass moments (to^-) and the combined mass- 
and-energy moments {n\) are almost identical. As de- 
scribed in the corresponding sections IV and V, the only 
differences regard the requirements needed to ensure a 
good resolution in the observables of interest. 

The analysis uses T(4S) — > BB events in which one 
of the B mesons decays to hadrons and is fully recon- 
structed (Bt a g), and the semileptonic decay of the recoil- 
ing B meson (-Brecon) is identified by the presence of an 
electron or muon. While this approach results in a low 
overall event selection efficiency of only a few per mille, it 
allows for the determination of momentum, charge, and 
flavor of the B mesons. 



A. Selection of Hadronic B-Meson Decays 

To obtain a large sample of -Btag-niesons, many exclu- 
sive hadronic decays B tag — > £)(*)y ± arc reconstructed 
[33]. The hadronic system consists of hadrons with 
a total charge of ±1. It is composed of n^^, nxK^, 



TABLE I: Summary of branching fractions of semilep- 
tonic decays B — > X C £~V used in MC simulations for 
neutral (B B o) and and charged (B B ±) B-meson decays. 
The values are taken from [21, 34-37]. Isospin symme- 
try is assumed to calculate the individual decay rates 
for B° and B mesons from their averaged measured 
branching fractions, taking into account the lifetime ratio 
1.071 ± 0.009 [21]. The sum of the exclusive decays is 
constrained to equal the total inclusive branching fractions 
for B -> X c i~V decays, B(B + -> XJ + u) = 10.89 ± 0.16 and 
B(B° -» XJ+u) = 10.15 ± 0.16 [34, 38]. 



Semileptonic Decay 


Bbo [%] 


B B ± [%J 


B 




Dl~V 


2.13 ± 0.14 


2.30 ± 0.16 


B 




DTv 


5.53 ± 0.25 


5.95 ± 0.24 


B 




Dxl~v 


0.50 ± 0.08 


0.54 ± 0.06 


B 




D* 2 t-V 


0.39 ± 0.07 


0.42 ± 0.08 


B 




D^rv 


0.43 ± 0.09 


0.45 ± 0.09 


B 


* 


D[£-V 


0.40 ± 0.20 


0.45 ± 0.20 


B 




D"tt£v 


0.40 ± 0.12 


0.20 ± 0.06 


B 




D ± tt£v 


0.19 ± 0.06 


0.40 ± 0.12 


B 




D*°n£v 


0.12 ± 0.04 


0.06 ± 0.02 


B 


* 


D* ± tt£v 


0.06 ± 0.04 


0.12 ± 0.04 



xi K aK s , and n w o7r° with n T + < 5, n^-o < 2, and 
n^o < 2, respectively. In total 1097 hadronic decay 
modes are reconstructed. 

The kinematic consistency of the B tag candidates is 
checked with two variables, the beam-energy-substituted 
mass 771es = \/s/4 — p§ and the energy difference AE = 
Eb — \/s/2. Here y/s is the total energy in the cm. frame, 
and ps and Eb denote the cm. momentum and cm. en- 
ergy of the -Btag candidate, respectively. The mass tties 
is measured with a resolution of 2.5 MeV/ c 2 , essentially 
independent of the Stag channel. We require AE = 
within three standard deviations, where one standard de- 
viation ranges between 10 and 30 MeV depending on the 
number of charged and neutral hadrons in the Btag can- 
didate. For each of the reconstructed hadronic modes the 
purity is estimated as the fraction of signal decays with 
toes > 5.27GeV/c 2 . We restrict the selection to hadronic 
modes with purities of at least 28% resulting in a selected 
Btag sample with an overal purity of 60%. On average we 
reconstruct B tag candidates with an efficiency of about 
0.4%. 



B. Selection of Semileptonic Decays 

Semileptonic decays are identified by the presence of 
one and only one electron or muon above a minimum mo- 
mentum p\ min measured in the rest frame of the B me- 
son. If not stated otherwise, p* t will denote in the follow- 
ing the lepton momentum measured in the -B-meson rest 
frame. Electrons are identified by combining information 
from the EMC, the DCH, and the DIRC. They are re- 
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quired to have a lab-frame momentum of p > 0.8 GcV/c 
and a polar angle in the range 0.41 < 9 < 2.54 rad. In 
this range, electrons are selected with 94% average effi- 
ciency and a hadron misidentification rate of the order of 
0.1%. Muon identification is mainly based on information 
obtained from the IFR. Muons are identified with an effi- 
ciency ranging between 60% for momenta p = 1 GeV/c in 
the laboratory frame and 75% for momenta p > 2 GcV/c. 
The misidentification rate ranges between 1% for kaons 
and protons and 3% for pions. Efficiencies and misidenti- 
fication rates are estimated from control samples of elec- 
trons, muons, pions, and kaons. We impose the condition 
QbQe < 0, where Qi is the charge of the lepton and Qb is 
the charge of the b quark of the B tag . This condition is 
fulfilled for primary leptons originating directly from the 
B decay, except for B°B° events in which flavor mixing 
has occurred. We require the total observed charge of the 
event to be |Q t ot| = l<3B tag + <2B rccoil | < 1, allowing for 
a charge imbalance in events with low momentum tracks 
or photon conversions. In cases where only one charged 
track is present in the reconstructed X c system, the total 
charge in the event is required to be zero. 



C. Reconstruction of the Hadronic System 

The hadronic system X c in the decay B — ► X c l~v 
is reconstructed from charged tracks and energy de- 
posits in the calorimeter that are not associated to the 
Stag or the charged lepton. We ignore tracks and en- 
ergy deposits in the calorimeter which are compatible 
with the hypothesis of being reconstruction artifacts, 
low-energy beam-generated photons or calorimeter de- 
posits originating from hadronic showers. Each track 
is assigned a specific particle type, either ft , K^, or 
7r ± , based on combined information from the different 
BABAR subdetectors. Few events containing single pro- 
tons are kept in the selection but removed later on in 
the background removal procedure. The four-momentum 
Px c of the reconstructed hadronic system is obtained 
from the four-momenta of the reconstructed tracks P^trk 
for the given mass assignment, and photons Pi :1 by 

p x c = Y*=i -fj.trk + Y^=i P hi- Tne hadronic mass is 
defined by m 2 x = P Xc . 

The four-momentum of the unmeasured neutrino is 
calculated from the missing four-momentum P m i ss = 
Pr(4S) ~ -Pfitag — Px c — Pt- Here, all four-momenta are 
measured in the laboratory frame. To ensure a well re- 
constructed hadronic system, we impose criteria on the 
missing energy, -E m i SS > 0.5 GeV, the missing momentum, 
Pmiss > 0.5GeV/c, and the difference of both quantities, 
l-Emiss - cp miss | < 0.5 GeV. 

We perform a kinematic fit exploiting the fact that 
B mesons are produced in a well-defined initial state 
e + e~ — > T(4S) — » BB, to improve the resolution and 
reduce the bias on the reconstructed quantities. The fit 
imposes four-momentum conservation, the equality of the 
masses of the two B mesons, and constrains the mass of 



the neutrino, -P miss = 0. The resulting (original) aver- 
age resolutions on the measurement of mx and n x are 
0.355 GeV/c 2 (0.425 GcV/c 2 ) and 1.05 GeV 2 (1.17 GeV 2 ), 
respectively. The average biases of mx and n 2 x after 
(before) the kinematic fit are found to be —0.096 GcV/c 2 
(-0.254 GeV/c 2 ) and -0.11 GeV 2 (-0.37 GeV 2 ), respec- 
tively. 

The background is composed of e + e~ — > qq(q = 
u,d,s,c) events (continuum background), T{AS) — > 
B + B~ or B°B° decays in which the B tag candidate is mis- 
takenly reconstructed from particles coming from both 
B mesons in the event (combinatorial background), and 
non-signal decays of the recoiling B r0C oii meson (residual 
background). Signal events where the hadronic system 
is not fully reconstructed are not considered as an addi- 
tional source of background. The effect of missing tracks 
and photons on the resolution of the kinematical quan- 
tities of interest is taken into account by applying the 
correction procedures described below. 

To quantify the amount of continuum and combinato- 
rial background in the toes signal region we perform a fit 
to the toes distribution of the -Btag candidates. We pa- 
rameterize the background using an empirical threshold 
function [39], 

« m ES Vl - x 2 e- x ( 1 -" 2 ), (2) 

dm E S 

where x = m ES /m ES ^ max , m E s,max = 5.289 GeV/c 2 is the 
kinematic endpoint approximated by the mean cm. en- 
ergy, and x is a free parameter defining the curvature of 
the function. The signal is parameterized with a modi- 
fied Gaussian function [40] peaked at the B-meson mass 
and corrected for radiation losses. The fit is performed 
separately for several bins in mx and n 2 x to account for 
changing background contributions. Figure 1 shows the 
toes distribution for p\ > 0.8GeV/c together with the 
fitted signal and background contributions. The shape 
of the continuum and combinatorial background as func- 
tion of mx and n 2 x is determined in a signal-free region 
of the toes sideband, 5.210 < toes < 5.255 GeV/c 2 . Its 
overall size in the toes signal region, toes > 5.27 GcV/c 2 , 
is determined by rescaling with the relative background 
contributions in the signal and sideband regions as de- 
termined by the fit. Signal and sideband region are sep- 
arated by 15 MeV/c 2 to avoid the leakage of signal events 
into the sideband region. 

Residual background is estimated from MC simula- 
tions. It is composed of charmless scmilcptonic decays 
B — > X u £~u, hadrons misidentified as leptons, secondary 
leptons from semileptonic decays of D(*> , Df mesons or r 
cither in B°B° mixed events or produced in b — ■> ccs tran- 
sitions, as well as leptons from decays of J/ip and ip(2S). 
The branching fractions of the individual simulated back- 
ground decays are scaled to agree with measurements 
[21, 34, 41, 42]. The overall simulated background spec- 
trum is normalized to the number of -Btag events in data. 
We verify the normalization and the shape using an inde- 
pendent data control sample with inverted lepton charge 
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FIG. 1: The toes spectrum of -Stag decays accompanied by a 
lepton with p| > 0.8GeV/c. The fit functions for the sum of 
signal and background (solid line) and the background (red 
dashed line) are overlaid. The crossed area shows the pre- 
dicted background under the Stag signal. The background 
control region in the toes sideband is indicated by the hatched 
area. 



correlation, QbQt > 0. 




IV. HADRONIC-MASS MOMENTS 

We present measurements of the moments (m x ), with 
k = 1, ... 6, of the hadronic-mass distribution in scmilcp- 
tonic £?-mcson decays B — > X c l~~D. The moments are 
measured as functions of the lower limit on the lepton 
momentum p* t min between 0.8GeV/c and 1.9GeV/c, cal- 
culated in the rest frame of the B meson. 



A. Selected Event Sample 

We find 19, 212 events with p* t > 0.8GcV/c, composed 
of 15, 085 ± 146 signal events above a combinatorial and 
continuum background of 2, 429 ± 43 events and residual 
background of 1, 696±19 events. Signal decays amount to 
79% of the selected event sample. For p* t > 1.9 GcV/c, wc 
find in total 2, 527 events composed of 2, 006 ± 53 signal 
events above a background of 271 ± 17 and 248 ± 7 com- 
binatorial/continuum and residual events, respectively. 
Figure 2 shows the mx distributions after the kinematic 
fit together with the extracted background shapes for 
p* e > 0.8GeV/c and p* e > 1.9GeV/c. 

B. Extraction of Moments 

To extract unbiased moments (m x ), we apply correc- 
tions to account for effects that distort the measured 
mx distribution. Contributing effects are the limited ac- 
ceptance and resolution of the BABAR detector resulting 
in unmeasured particles and in misreconstructed ener- 
gies and momenta of particles. In addition, there are 
contributions from measured particles not belonging to 



FIG. 2: Hadronic-mass spectra after the kinematic fit for lep- 
ton momenta p\ > 0.8GeV/c (top) and p\ > 1.9GeV/c (bot- 
tom) together with distributions of combinatorial background 
and background from non-SB decays (red, hatched area) as 
well as residual background (blue, crossed area). The two 
background histograms are plotted on top of each other. 



the hadronic system, especially photons originating from 
FSR of the primary leptons. These photons are included 
in the measured X c system and thus lead to a modified 
value of its mass; they also lower the momentum of the 
primary lepton. Both effects are included in our correc- 
tion procedure. 

We correct the kinematically-fitted value of m\ of each 
event by applying correction factors on an event-by-event 
basis using the observed linear relationship between the 
moments of the measured mass (m x reco ) and the mo- 
ments of the true mass (m x t ) in MC spectra. The 
correction factors are determined from MC simulations 
by calculating moments (m x reco ) and (m x true ) in sev- 
eral bins of the true mass mxjrue and fitting the ob- 
served dependence with a linear function, referred to as 
calibration function in the following. 

We find that the bias of the measured moments 
(m x reco ) is not constant over the whole phase space. 
Therefore, we derive the calibration functions in three 
bins of the particle multiplicity Nx c in the X c system, 
three bins of -Emiss — cp m i ss , as well as in twelve bins of 
p* t , each with a width of lOOMeV/c. Due to the limited 
number of generated MC events, the binning in Nx c and 
Emiss - cp miss is not used for p* i min > 1.7GeV/c. Over- 
all we construct 84 calibration functions for each order 
of moments. The obtained calibration functions allow a 
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FIG. 3: Examples of calibration functions for {m 2 x ) in bins of Nx a , -E m i ss — cp m j aa and p\ . Shown are the extracted moments 
i m x,reco) versus the true moments {m x true ) for 0.8 < p| < 0.9GeV/c (•), 1.4 < p\ < 1.5GeV/c (o), and p\ > 1.9GeV/c (■). 
The results of fits of linear functions are overlaid as solid lines. A reference line with {m 2 x jreco ) = (m x true ) is superimposed 
(dashed line). The calibration function for p\ > 1.9 GeV/c is constructed independent of Nx c and Emi BB — cp m i ss . It is plotted 
in each of the bins for comparison only. 



consistent extraction of moments for events containing 
an electron or a muon. Figure 3 shows examples of cali- 
bration functions for the moment {m x } in three bins of 
p\ as well as in nine bins of -Emiss — cp m i SS and N Xc ■ 

For each data event i, the corrected mass m x calib i is 
calculated by inverting the linear function, 



k _ m X,reco,i - ^(^miss ~ CPmiss, N Xc , k,pf) 

.calib. i t-w TTi aT t *\ ' 

B{E miss - qjnuss) N Xc ,k,Pi) 

(3) 

where A is the offset and B is the slope of the calibra- 
tion function. Background contributions are removed by 
applying a weight factor Wi to each corrected hadronic 
mass m x calib - , where the weight is the expected frac- 
tion of signal events in the corresponding region of the 
mx,reco spectrum in Fig. 2. The expression used to cal- 



culate the moments is the following: 

E Wi{m x )m k X calib l 
(m x ) = — — x C ca i(p* e , k) X Ctrueipl, k), 

i 

(4) 

with N ev the total number of selected events. The fac- 
tors Ccai and Ctme depend on the order k and the mini- 
mum lepton momentum p* t ■ of the measured moment. 
They are determined in MC simulations and correct for 
the residual small biases observed after the calibration. 
The factors C ca i account for the bias of the applied cor- 
rection method and are typically ranging between 1.01 
and 1.06 for k = 1 ... 5. Larger bias corrections C ca i are 
observed for (m x ) ranging between the extremes 0.902 
and 1.109. The residual bias-correction factor Ctme ac- 
counts for differences in selection efficiencies for different 
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hadronic final states and FSR that is included in the mea- 
sured hadron mass and distorts the measurement of the 
lepton's momentum. Our correction procedure results in 
moments which are within systematic uncertainties free 
of photon radiation. The correction Ct rue is estimated in 
MC simulations and typically ranges between 0.994 and 
1.007. For the moments (m 5 x ) and (m x ), slightly higher 
correction factors are determined, ranging between 0.990 
and 1.014 for (m x ) and 0.986 and 1.024 for (m x ). 

This correction procedure is verified on a MC sample 
by applying the calibration to measured hadron masses of 
individual scmileptonic decays, B — > D£~Ts, B — > D*£~17, 
four resonant decays B — > D**£V, and two nonresonant 
decays B -»■ D^iriU. Fi gure 4 shows the corrected 
moments (m 2 x ) and (m x ) as functions of the true mo- 
ments for minimum lepton momenta p\ > 0.8GeV/c. The 
dashed line corresponds to (m x calib ) = (m xtrue ). The 
calibration reproduces the true moments over the full 
mass range. 



C. Systematic Uncertainties and Tests 

The main systematic uncertainties are associated with 
the modeling of hadronic final states in scmileptonic B- 
meson decays, the bias of the calibration method, the 
determination of residual background contributions, the 
modeling of track and photon selection efficiencies, and 
the identification of particles. The total systematic un- 
certainty is estimated by adding in quadrature five con- 
tributions, as described below. Tables A.I and A. II list 
the individual contributions to the systematic errors of 



the measured moments 



x < 



with k = 1 ... 6 for 



mum lepton momenta ranging from 0.8 to 1.9GeV/c. 



1. MC Statistics 

The effect of limited MC statistics on the extracted 
moments is evaluated using parameterized MC experi- 
ments. To study the effect on the calibration curves, 
the parameters of the fitted first-order polynomials are 
randomly varied within their uncertainties including cor- 
relations and new sets of moments are extracted. The 
overall uncertainty is determined by repeating this pro- 
cedure 250 times and taking the r.m.s. of the distribution 
of the moments as the systematic uncertainty. 

To estimate the effect of limited MC statistics in the 
residual background determination a similar method is 
applied by varying the parameters of the fit to the /ties 
distribution within their errors including correlations. 



2. Simulation- Related Effects 

We correct for differences between data and MC sim- 
ulation in the selection efficiencies of charged tracks and 
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FIG. 4: Calibrated (•) and uncorrected (□) moments (m 2 x ) 
(left) and (rn\) (right) of individual hadronic modes for 
lepton momenta p} > 0.8 GeV/c. A reference line with 
{■m x ,caiib} = {mx.true) is superimposed. 



photons, as well as identification efficiencies and misiden- 
tification rates of various particle types. The corrections 
are extracted from data and MC control samples. 

The systematic uncertainties of the photon selection 
and track finding efficiencies are determined studying in- 
dependent control samples. Their impact on the mea- 
sured moments has been evaluated by randomly exclud- 
ing neutral or charged candidates from the X c system 
with probabilities corresponding to the uncertainties of 
the efficiency extraction methods. The uncertainty of the 
photon selection efficiencies is found to be 1.8% per pho- 
ton independent of energy, polar angle, and multiplicity. 
The uncertainty in track finding efficiencies consists of 
two parts. For each track, we add in quadrature 0.8% 
systematic uncertainty and the statistical uncertainty of 
the control samples that depend on energy and polar an- 
gle of the track as well as the multiplicity of tracks in the 
reconstructed event. 

The systematic uncertainty on the misidentification of 
7T^~ mesons as leptons is found to affect the overall nor- 
malization of the corresponding background spectra by 
8%. The influence on the measured moments is estimated 
by varying the corresponding background within its un- 
certainty. The observed variation of moments is taken as 
a systematic uncertainty. 

The impact of mismodeling FSR simulated with 
PHOTOS [26] is estimated by calculating moments from 
data using a set of calibration curves constructed from 
events simulated without FSR photons. The theoretical 
uncertainty associated with the calculations included in 
PHOTOS is conservativley assumed to be of the order of 
20%. The systematic uncertainty connected to the mis- 
modeling of FSR photons is therefore estimated to be 
20% of the observed difference between the nominal mo- 
ments and those from the MC simulation without FSR 
photons. 

A significant fraction of the low-energy photons de- 
tected in the calorimeter are beam related. We check 
the impact of low-energy photons by removing EMC sig- 
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nals with energies below 100 McV from the reconstructed 
hadronic system. The effect on the measured moments 
is found to be negligible. 

The stability of the result under variation of the se- 
lection criteria on E m i SS — qj m iss is tested by varying 
the applied cut between |E m iss — cpmissl < 0.2 GeV and 
l-E-miss — cpmissl < 1.4 GeV. For all measured moments, 
the observed variation is well covered by other known sys- 
tematic detector and MC simulation effects. Therefore, 
no contribution is added to the systematic uncertainty. 

3. Extraction Method 

The systematic uncertainty of the calibration bias cor- 
rection C ca i is estimated to be {C ca i — l)/2- 

4- Background Determination 

The branching fractions of background decays in the 
MC simulation are scaled to agree with the current mea- 
surements [21, 34, 41, 42]. The associated systematic 
uncertainty is estimated by varying these branching frac- 
tions within their uncertainties. At low p| min , most of 
the studied background channels contribute to the sys- 
tematic uncertainty equally, while at high p\ min , the sys- 
tematic uncertainty is dominated by background from 
decays B — > X u i~v. Contributions from J/ip and ip(2S) 
decays are found to be negligible. 

The uncertainty in the combinatorial -B tag background 
determination is estimated by varying the lower and up- 
per limits of the sideband region in the mES distribution 
up and down by 2.5MeV/c 2 . The observed effect on all 
hadronic-mass moments is found to be negligible. 

5. Modeling of Signal Decays 

The uncertainty of the calibration method with re- 
spect to the chosen signal model is estimated by changing 
the composition of the simulated inclusive hadronic spec- 
trum. The dependence on the simulation of high mass 
hadronic final states is estimated by constructing calibra- 
tion functions only from MC simulated hadronic events 
with hadronic masses mx,true < 2.5GeV/c 2 , thereby re- 
moving the high mass tail of the simulated hadronic- 
mass spectrum. The model dependence of the calibration 
method is found to be a small contribution to the total 
systematic uncertainty. 

We estimate the model dependence of the residual bias 
correction Ct rU e by changing the composition of the inclu- 
sive hadronic spectrum, i.e. omitting one or more decay 
modes. 

We study the effect of differences between data and 
MC simulation in the multiplicity and E m iss — cp m iss dis- 
tributions on the calibration method by changing the bin- 
ning of the calibration functions. The observed variation 



of the results are found to be covered by the statistical 
uncertainties of the calibration functions, and no contri- 
bution is added to the total systematic uncertainty. 



6. Stability of the Results 

The stability of the results is tested by dividing the 
data into several independent subsamplcs: B^ and E°, 
decays to electrons and muons, different run periods 
of roughly equal sample sizes, and two regions in the 
E miss - cp miss spectrum, -0.5 < E miss - cp miss < GeV 
and < -E m iss — cp m ; ss < 0.5 GeV, characterized by dif- 
ferent resolutions of the reconstructed hadronic system. 
No significant variations are observed. 



D. Results 

The measured hadronic-mass moments (m x ) after ra- 
diative correction with k = 1 ... 6 as functions of the 
minimum lepton momentum p* e min are shown in Fig. 5. 
All measurements are correlated since they share subsets 
of selected events. Tables A.I and A. II summarize the 
numerical results. In most cases wc find systematic un- 
certainties that exceed the statistical uncertainty by a 
factor of 2.5. The correlation matrix for the moments is 
given in the EPAPS document [43]. 



V. MOMENTS OF THE COMBINED 
MASS-AND-ENERGY SPECTRUM 

The measurement of moments of the observable n 2 x , a 
combination of the mass and energy of the inclusive X c 
system, as defined in Eq. (1), is theoretically motivated 
and is expected to allow a more reliable extraction of the 
higher order HQE parameters fi 2 and (? D [20]. 

We present measurements of the moments (n x ), (n x ), 
and (n x ) for different minimum lepton momenta between 
0.8GeV/c and 1.9GeV/c in the E-meson rest frame. 



A. Event Selection 

Due to the structure of the variable difference 
of two measured values, its measured resolution and bias 
are worse than for the mass moments. Also, the sensitiv- 
ity to cuts on Emiss — cpmiss increases. The average res- 
olution of n 2 x after the kinematic fit for lepton momenta 
greater than 0.8 GeV/c is measured to be 1.05 GeV 2 with 
a bias of -0.11 GeV 2 . We therefore introduce stronger 
requirements on the reconstruction quality of the event. 
We tighten the criteria on the neutrino observables by 
requiring Emiss — cpmiss to be between —0.2 and 0.3 GeV. 
Due to the stronger requirement, the individual variables 
Emiss and p m i ss have less influence on the resolution of 
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FIG. 5: Radiation-corrected hadronic-mass moments (m^-) with k = 1...6 for different selection criteria on the minimum 
lepton momentum Pf jm i n . The inner error bars correspond to the statistical uncertainties while the full error bars correspond 
to the total uncertainties. The moments, as well as their values for different Pe ilrllll , are highly correlated. 



the reconstructed hadronic system. Therefore, the re- 
quirements on the missing energy and the missing mo- 
mentum in the event are relaxed to E m i SS > GeV and 
Pmiss > GeV/c, respectively, as these requirements do 
not yield significant improvement on the resolution of 
n x , and do not increase the ratio of signal to background 
events. 

For p* ( > 0.8GcV/c and 1.9GeV/c, there remain 
10, 053 ± 142 and 1, 626 ± 52 signal events, respectively 
Background events make up 22 % of the final event sam- 
ple with p\ > 0.8 GeV/c. The background is composed of 
12 % continuum and combinatorial background and 10 % 
decays of the signal B meson other than the semilep- 
tonic decay B — > X c l~v. Combinatorial and contin- 
uum background is removed using the sideband of the 
toes distribution, as described in section IIIC. The 
residual background events, containing a correctly recon- 
structed Stag meson, are removed using MC simulations. 
The dominant sources are pions misidentified as muons, 
B — > X u £~v decays, and secondary semileptonic decays 
of D and D s mesons. 

The measured r? x spectra for p| min = 0.8GeV/c and 
P*i min = l-9GeV/c are shown together with the back- 
ground distributions in Fig. 6. 



B. Extraction of Moments 

The extraction of unbiased moments (n x ) from the 
measured n x spectra follows a calibration procedure sim- 



ilar to the one used to extract the hadronic-mass mo- 
ments as described in Section IV B. The linear calibra- 
tion functions 

k _ n X,y eCO ~ ^(ffmiss ^ Cp miss , N Xc ,k,p* t ) 

B{E miss - cp miss , N Xc ,k,p e ) 

for k = 2, 4, 6 are derived from MC samples in three bins 
of -E m iss — cp m i ss and three bins of the ^-system multi- 
plicity Nx c for each of the 12 lepton momentum bins of 
lOOMeV/c width. Because of differences in events con- 
taining electrons and muons, we also derive separate cali- 
bration functions for these two classes of events. Overall, 
we determine 216 linear calibration functions. The cali- 
bration again includes the effects of FSR photons which 
not only modify mx and p* t , but also Ex ■ 

We have verified that applying the calibration proce- 
dure on MC samples of individual exclusive B — > X c l~v 
modes allows to reproduce the generated moments, as 
shown in Fig. 7. Small biases remaining after calibration 
are of the order of 1 % for (n x ) and of few percent for 
(n x ) and (n%). 

Background contributions are removed by applying 
rt^-dependent weight factors Wi(n x ) on an event-by- 
event basis, leading to the following expression for the 
determination of the moments: 

E Wi{n x ) n Xcahhi 

(4) = ^ x c ^^- ( 6 ) 

E Wi{n x ) 

i=l 
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FIG. 6: Spectra of n x after the kinematic fit together with 
distributions of combinatorial background and background 
from non-SB decays (red, hatched area) as well as residual 
background (blue, crossed area) for different minimum lepton 
momenta (a) p? jm i n = 0.8GeV/c and (b) pl iinin = 1.9GeV/c. 
The two background histograms are plotted on top of each 
other. 



a 5 

V 

1 3 
V 2 
1 

0, 



"slope = 


1.002 ±0.009 


* 


^offset = 


-0.015+0.013 


/ 

/ 




* 

/'to 










: D* 






TJ P 

: 9 

\,. . . . i . . 




(a) 




. . i . . . . i . . . . i . . . 





3 



> 

02 



slope = 0.974 ±0.021 
5foff'sct = -0.020 ± 0.026 J 



] 2 
.1.5 
1 

0.5 



D/ 
* 

□ 



D ' 

if 



' 



(b) 



1 



2 3 



4 5 6 
ue > [GeV 2 ] 



0.5 1 



1.5 

2 

X.true 



<n 



2 2.5 3 
> [GeV 2 ] 



FIG. 7: Example of the calibration verification procedure for 
different minimum lepton momenta (a) p! m j n = 0.8GeV/c 
and (b) p} m i n = 1.7GeV/c. Moments {n 2 x ) of exclusive modes 
on simulated events before (□) and after (•) calibration are 
plotted against the true moments for each mode. The dotted 
line shows the result of a fit to the calibrated moments, the 
resulting parameters are given. 



The bias correction factors C(p^,k), depending on the 
minimum lepton momentum and the order of the ex- 
tracted moments, are determined by MC simulations; 
they combine the two factors C ca i and C trU e &s described 
in Section IV B. 



C. Systematic Uncertainties and Tests 



We consider the same five sources of systematic uncer- 
tainties as for the mass moments described in Sections 
IV CI to IV C 5: MC statistics, simulation-related ef- 
fects, extraction method, background determination, and 
modeling of signal decays. The individual contributions 
to the systematic error, listed in Table A. Ill, are esti- 
mated following procedures essentially identical to those 
described for the mass moments. 

Because of the tighter cut on -E m i ss — cp m ; ss , the sys- 
tematic uncertainty associated with this criterion is esti- 
mated in a different way. Wc first keep the lower limit 
fixed to the nominal value and vary the upper limit to 
0.3GeV/c to 0.25GcV/c, 0.4GeV/c, and 0.5GcV/c. Then 
we fix the upper limit to its nominal value and vary the 
lower limit to — 0.3GeV/c and — O.lGcV/c. The mean 
of the observed differences in the measured moments on 
data is taken as systematic uncertainty. 

In the third study, we include the uncertainty from the 
binning of the calibration function in the multiplicity of 
the X c -system. For the choice of the calibration function, 
we randomly increase the measured multiplicity of the X c 
system by one with a probability of 5% corresponding 
to the observed difference between MC and data. The 
uncertainty in the bias-correction factor C(pl,k) is con- 
servatively estimated as half of the applied correction. 

Varying the branching fractions of the exclusive signal 
modes in the MC simulation has, in agreement with the 
mass-moment studies, a very small impact on the mea- 
sured combined moments. Also, no significant variations 
of the results are observed when splitting the data sample 
into the same subsamples as for the mass moments. 



D. Results 



Figure 8 shows the results for the moments (n^-), (nj^), 
and {n\) as a function of the minimum lepton momen- 
tum p\ min . The moments are highly correlated due to the 
overlapping data samples. The full numerical results and 
the statistical and the estimated systematic uncertain- 
ties are given in Table A. III. The systematic covariance 
matrix for the moments of different order and with differ- 
ent cuts on p\ min is built using statistical correlations. 
This correlation matrix for the moments is given in the 
EPAPS document [43]. 

A clear dependence on the minimum lepton momen- 
tum is observed for all moments, due to the increasing 
contributions from higher-mass final states with decreas- 
ing lepton momentum. In most cases we obtain system- 
atic uncertainties slightly exceeding the statistical uncer- 
tainty. 
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FIG. 8: Radiation-corrected moments (a) (n%), (b) (n%), and (c) (n%) for different cuts on the lepton momentum p|. The 
inner error bars correspond to the statistical uncertainties while the full error bars correspond to the total uncertainties. The 
moments are highly correlated. 



VI. MOMENTS OF THE ELECTRON-ENERGY 
SPECTRUM 

Moments of the electron-energy spectrum for semilcp- 
tonic decays B — > X c e~v averaged over charged and 
neutral B mesons have been measured in a data sam- 
ple of 51 x 10 6 r(AS) BB decays [9]. In the following, 
we present an overview of this analysis and update the 
results by using more recent measurements [21, 41] of 
branching fractions of background processes. 

In multi-hadron events as defined in [9], BB events 
are selected by requiring a semilcptonic B decay with an 
identified electron (etag), with charge Q(e t ag) and a mo- 
mentum 1.4 < p* < 2.3GcV/c, measured in the T(4S) 
rest frame. These events constitute a tagging sample 
used as normalization for the branching fraction. A sec- 
ond electron e S i g , for which we require p* > 0.5GeV/c, 
is assigned either to the unlike-sign sample if the tagged 
sample contains an electron with Q(etag) = — <3(e s i g ) or 
to the like-sign sample if <5(e ta g) = Q(e s i g )- In events 
without B"B" mixing, primary electrons from semilcp- 
tonic B decays belong to the unlike-sign sample while 
secondary electrons contribute to the like-sign sample. 
Secondary electrons originating from the same B as the 
e tag are removed from the unlike-sign sample by the re- 
quirement 

cosa* > 1.0-p*c/GeV and cos a* > -0.2, (7) 

where a* is the angle between the two electrons in 
the T(4S) rest frame. Corrections for the small resid- 
ual background of unlike-sign pairs originating from the 
same B fulfilling this requirement arc described in [9]. 
Additional background corrections for electrons from 



j/4> 



decays, continuum events, photon con- 



versions, 7r — > e + e~7 Dalitz decays, and misidentificd 
hadrons are also described in [9], Figure 9 shows the 
electron-momentum spectra together with the contribu- 
tions of the backgrounds. 

Further backgrounds arise from decays of r leptons, 
charmed mesons produced in b — > ccs decays, and J/tp 



or tp(2S) — > e + e~ decays with only one detected elec- 
tron. We also need to correct for cases where the tagged 
electron does not originate from a semilcptonic B de- 
cay. These backgrounds are irreducible. Their contribu- 
tions to the three samples - single electrons, like-sign, and 
unlike-sign pairs - are estimated from MC simulations, 
using the ISGW2 model [31] to describe semilcptonic D 
and Z? s -meson decays. As an important update to the re- 
sults in [9] , the branching fractions of these backgrounds 
are recalculated to match the recent measurements [211. 
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FIG. 9: Previously measured momentum spectrum (points) 
[9] and estimated backgrounds (histograms) for electron can- 
didates in (a) the unlike-sign sample, and (b) the like-sign 
sample. The background spectra are updated wrt. the previ- 
ous publication with more recent branching-fraction measure- 
ments [21, 41]. 
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As in [9], we calculate B(D S -> Xev) = (7.79 ± 0.19)% 
from B(D° — * Xev) and B(D + — ► Xev), assuming 
r(L> s -» Xev) = T(D -> Xei/). Using -> 
£>+X) = (8.3 ± 0.8)% [21] the branching fraction of 
B°' + — > £>+ — > e + decays, where the D s originates 
from fragmentation of the W boson, is (0.65 ± 0.06)%. 
Using the inclusive branching fraction measurement of 
— ► D°' + X decays reported in [41], we arrive at 
B(B°>+ -> L>°-+ -> e+) = (0.93 ± 0.11)%. To estimate 
the contribution of electrons from r decays, we consider 



the cascades B — > r — > e and B — > Z) s — > r — > e, with 
branching fractions taken from [21]. The rates for the 
decays B — > <//(/> — ► e + e~ and S — > i/)(2S) — > e + e~ are 
also adjusted to [21]. 

After the like- and unlike-sign samples have been cor- 
rected for electron identification efficiency, these irre- 
ducible background spectra are subtracted. To account 
for B°B° mixing, we determine the number of primary 
electrons in the i-th. p*-bin from the like-sign and unlike- 
sign pairs as 



J 



(l-/oXo)-(l-p)(l-/o) 



(l-2/oxo)-(l-p)(l-/o)(l-/oXo) 



Nl 



foXo 



(l-2/o X o)-(l-p)(l-/o)(l-/oXo) 



(8) 



r 



where xo = 0.1878 ± 0.0024 [21] is the B°B° mixing pa- 
rameter, f = B{T(AS) -> B°B°) = 0.491 ± 0.007 [21], 
and p = B{B+ -► D° -> e -)/B(B° D~ e~) = 
(0.744 ±0.06) [21]. The parameter e l a , is the efficiency of 
the additional requirement for the unlike-sign sample as 
defined in Eq. (7). The spectrum obtained from Eq. (8) 
is corrected for the effects of brcmsstrahlung in the de- 
tector material using MC simulation. Figure 10 shows 
the resulting spectrum of primary electrons. 

We determine the partial branching fraction as 
(Hi N b^c,u) / ( N tng e evt (-cuts), where i runs over all bins 
with E e > Eq. For the background-corrected number 
A tag of tag electrons we find A tag = (3617±4±22) x 10 3 , 
where the uncertainties are statistical and systematic, re- 
spectively. The parameter e evt — (98.9 ± 0.5)% refers 
to the relative efficiency for selecting two-electron events 
compared to events with a single et ag , and e cu t s = 
(82.8 ± 0.3)% is the acceptance for the signal electron 
for E = 0.6 GeV. The result is 

B(B -> Xev{j), E e > 0.6 GeV) 

= (10.30 ±0.06 ±0.21)%, 

where the errors correspond to the statistical and sys- 
tematic uncertainties, respectively. 

In the B-meson rest frame, we define Ri(Eo,n) as 
f™(E e - p) l (dT/dE e )dE e , and measure the first mo- 
ment Mi (S ) = Ri(E o ,0)/R a (E o ,0), the central mo- 
ments M n (E ) = R n (E , M 1 (E ))/R (Eo,0) for n = 2, 3 
and the partial branching fraction B(E n ) = tb Ro(E ,0), 
where tb is the average lifetime of charged and neutral 
B mesons. The calculation of the moments is done as in 
[9] and includes corrections for charmless semilcptonic de- 
cays, the movement of the B mesons in the cm. frame, 
biases due to the event selection criteria, and binning 
effects. The spectra and moments presented are those 
of B — > X c ev(-y) decays with any number of photons. 
Since current theoretical predictions on the lepton-encrgy 
moments do not incorporate photon emission, we also 
present a second set of of moments with corrections for 



the impact of QED radiation using the PHOTOS code [26] . 

Figure 11 shows the moments of B — + X c ev(j) decays 
as a function of Eo, and Table II lists the main systematic 
errors for Eq = 0.6 and 1.5 GeV. The complete listing of 
all moments and the full correlation matrix, with and 
without PHOTOS corrections can be found in [43]. 



VII. DETERMINATION OF \V cb \ AND THE 
QUARK MASSES m b AND m c 

At the parton level, the weak decay rate for b — * civ 
can be calculated accurately; it is proportional to \V c b\ 2 
and depends on the quark masses nib and m c . To relate 
measurements of the semileptonic B-meson decay rate to 
\V c b\, the parton-levcl calculations have to be corrected 
for effects of strong interactions. Heavy Quark Expan- 
sions (HQEs) [44-46] have become a successful tool for 
calculating perturbative and nonperturbative QCD cor- 
rections [47-51] and for estimating their uncertainties. 
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FIG. 10: Electron-momentum spectrum from B — ► Xev(-y) 
decays in the T(4S) frame after correction for efficiencies and 
bremsstrahlung in the detector, with combined statistical and 
systematic errors. 
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FIG. 11: Measured branching fraction (a) and moments Mi (b), M2 (c), and M3 (d) of the inclusive electron-energy spectrum 
of B —> Xcevi^y) decays as a function of the cutoff energy Eq in the B-meson rest frame. 



In the kinetic-mass scheme [11, 20, 52-55], these 
expansions in 1/rrib and the strong coupling constant 
ctsijnb) to order C(l/m^) contain six parameters: the 
running kinetic masses of the b and c quarks, m(,(/j) and 
m c (/j), and four nonperturbative parameters. The pa- 
rameter n denotes the Wilson factorization scale that 
separates effects from long- and short-distance dynam- 
ics. The calculations are performed for \i = 1 GeV 
[56]. It has been shown that the expressions for the 
moments have only a small scale dependence [17]. We 
determine these six parameters and | V c b | from fits to mo- 
ments of the hadronic-mass, combined mass-and-encrgy, 
and electron-energy distributions in semileptonic B de- 
cays B — > X c l~~D and moments of the photon-energy 
spectrum in decays B — > X s j [14-16]. 

In the kinetic-mass scheme the HQE to C(l/m^) for 
the rate T$l of semileptonic decays B — > X c £~v can be 
expressed as [11] 



-£-±\V cb \ 2 (l + A ew )A pert (r,p) 



Mr) 1 



,,2 ,,2 , Pd+pI 



(9) 



- 2(1 -r) 
+ 0(l/mt) 



4 A*g 



Pd+Pl 



cm 



C 777.1, 



The leading nonperturbative effects arise at 0(1 /to 2 ) 
and are parameterized by /i 2 (/j) and /jg (/j) , the expecta- 
tion values of the kinetic and chromomagnctic dimension- 
five operators. At 0(l/mf), two additional parameters 
enter, p 3 D {p) and p\g{p)i the expectation values of the 
Darwin and spin-orbit dimension-six operators, respec- 
tively. The ratio r = m 2 /m 2 enters in the tree level 



phase space factor zo(r) = 1 — 8r + 8r 3 — r 4 — 12r 2 lnr 
and in the function d(r) = 81nr + 34/3 - 32r/3 - 8r 2 + 
32r 3 /3 — 10r 4 /3. The factor 1 + A ew accounts for elec- 
troweak corrections. It is estimated to be 1 + A ew = 



(1 + a/ir In M z /m b f 



1.014, where a is the electro- 



magnetic coupling constant. The quantity A pert accounts 
for perturbative contributions and is estimated to be 
A pert {r, fj.) ^ 0.908 [11]. 

The performed fit uses a linearized expression for the 
dependence of \V c b\ on the values of heavy-quark param- 
eters, expanded around a priori estimates of these pa- 
rameters [111: 



Ib{b -> xj-v) 1.55 



(10) 



0.0417 V °- 1032 TB 

x[l + 0.30(a s (m&)-0.22)] 
x [1 - 0.66(to 6 - 4.60) + 0.39(m c - 1.15) 
+0.01304 - 0.40) + 0.09(pf> - 0.20) 
+0.05(/4 - 0.35) - 0.01(p| s + 0.15)]. 

Here m& and m c are in GeV/c 2 and all other param- 
eters of the expansion are in GeV fc ; tb refers to the 
average lifetime of B mesons produced at the T(4S), 
measured in picoseconds. HQEs in terms of the same 
heavy-quark parameters are available for hadronic-mass, 
combined mass-and-energy, electron-energy, and photon- 
energy moments. Predictions for those moments arc ob- 
tained from an analytical calculation [57]. We use these 
calculations to determine \V cb \, the total semileptonic 
branching fraction B(B — > X c £~v), the quark masses mt, 
and m c , as well as the heavy-quark parameters /j 2 , /j^., 
and p\ s , from a simultaneous x 2 fit to the measured 
moments and partial branching fractions, all as functions 
of the minimum lepton momentum p* t min and minimum 
photon energy -E 7 , min . 
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TABLE II: Results and breakdown of the systematic errors for the partial branching fraction B — tb J^' g (dT/dEe) dE e , and the 
moments Mi, M2, and M3 for B — > X c ev and B — » X c ev(^) in the B-meson rest frame for two values of Eq. Changes wrt. the 
previously published results [9] are due to updated background branching fractions [21, 41] and indicated by (+). 





6[io^] 


Mi[ 


MeV] 


M 2 [10~ 


" a GeV^] 


Af 3 [10- a 


GeV 11 ] 


£ [GeV] 


0.6 


1.5 


0.6 


1.5 


0.6 


1.5 


0.6 


1.5 


Breakdown of systematic errors 


Conversion and Dalitz pairs 0.028 


0.001 


1.5 


0.02 


0.6 


0.00 


0.04 


0.00 


e identification efficiency 


0.150 


0.044 


2.5 


0.30 


0.6 


0.07 


0.27 


0.08 


e from same B 


0.019 


0.000 


1.3 


0.00 


0.6 


0.00 


0.05 


0.00 


B -► D s -► e (f) 


0.024 


0.000 


1.4 


0.01 


0.5 


0.00 


0.03 


0.00 


B - D - e (f) 


0.035 


0.000 


2.2 


0.00 


1.0 


0.00 


0.03 


0.00 


B - r -* e (f) 


0.027 


0.001 


1.2 


0.04 


0.3 


0.00 


0.10 


0.00 


e from J/V> or ^(25*) (f) 


0.002 


0.001 


0.0 


0.01 


0.0 


0.01 


0.00 


0.00 


Secondary tags 


0.052 


0.011 


1.6 


0.06 


0.6 


0.00 


0.05 


0.00 


X 


0.022 


0.011 


0.9 


0.01 


0.3 


0.00 


0.03 


0.00 


Tracking efficiency 


0.083 


0.033 


1.0 


0.06 


0.3 


0.02 


0.06 


0.00 


Bremsstrahlung correction 


0.011 


0.028 


1.9 


0.43 


0.0 


0.05 


0.18 


0.00 


Event selection 


0.052 


0.024 


0.6 


0.14 


0.0 


0.03 


0.07 


0.01 


b — > u subtraction (f) 


0.031 


0.020 


0.8 


0.83 


0.4 


0.32 


0.14 


0.12 


B momentum correction 


0.000 


0.005 


0.0 


0.19 


0.1 


0.10 


0.04 


0.02 


Ntag normalization 


0.068 


0.030 














Results 


Results for B — ► X ev(^j) 


10.08 


4.53 


1418.8 


1768.7 


146.1 


29.6 


-10.08 


2.04 


±(stat.) 


0.06 


0.03 


3.8 


1.9 


2.0 


0.8 


0.81 


0.44 


±(syst.) 


0.21 


0.08 


5.4 


1.0 


1.9 


0.4 


0.40 


0.15 


Results for B — > X c ei^ 


10.20 


4.78 


1437.6 


1773.7 


145.4 


30.1 


-12.04 


2.04 


±(stat.) 


0.06 


0.03 


4.0 


1.9 


2.3 


0.9 


0.91 


0.47 


±(syst.) 


0.22 


0.08 


5.7 


1.1 


2.1 


0.4 


0.40 


0.15 



Extraction Formalism 



parameter in the fit. It is given by 



The fit method designed to extract the HQE parame- 
ters from the measured moments has been reported pre- 
viously [17, 58]. It is based on a \ 2 minimization, 



M fi . 



M, 



HQE 



) T c to j( 



M, 



M 



HQE 



(11) 



The vectors M exp and Mhqe contain the measured mo- 
ments and the corresponding moments calculated by the- 
ory, respectively. Furthermore, the expression in Eq. (11) 



contains the total covariancc matrix C*. 



C c 



C 



HQE 



defined as the sum of the experimental C cxp and theoret- 
ical Chqe covariance matrices (see Section VII C). 

The total semileptonic branching fraction B(B — > 
X C £~T>) is extracted in the fit by extrapolating the mea- 
sured partial branching fractions B p * . (B — * X c £~~p) 
with > p\ min to the full lepton energy spectrum. Us- 
ing HQE predictions of the relative decay fraction 



R 



J P?,mi„ d Pe 1 



Pf,»i„ 



(12) 



JO dp* a Pi 

the total branching fraction can be introduced as a free 



B(B -> X c rv) = Pl ^\ -. (13) 

Using Eqs. (10) and (11) together with the measured av- 
erage B-meson lifetime tb and the total branching frac- 
tion, allows the calculation of \V c b\- 



|u cfc | 2 ocr SL = & 



xj- 



TB 



(14) 



Thereby, \V c b\ is introduced as an additional free param- 
eter to the fit. To propagate the uncertainty on tb prop- 
erly into the extracted result for \V c b\, tb is added as an 
additional measurement to the vectors of measured and 
predicted quantities, M exp and A^hqe- 

The nonperturbative parameters [£q and p\ s have 
been estimated from the B-B* mass splitting and heavy- 
quark sum rules to be p? G = (0.35 ± 0.07) GeV 2 and 
Pls = (-0.15 ±0.10) GeV 3 [17], respectively. Both pa- 
rameters are restricted in the fit by imposing Gaussian 
error constraints. 



B. Experimental Input 

The combined fit is performed on a subset of available 
moment measurements with correlations below 95% to 
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ensure the invertibility of the covariance matrix. Since 
the omitted measurements are characterized by high cor- 
relations to other measurements considered in the fit, 
they do not contribute significant additional information, 
and the overall sensitivity of the results is not affected. 
Choosing a different subset of moments gives consistent 
results. We perform two fits to the following set of mea- 
sured moments, thereby including either the hadronic- 
mass moments or the moments of the combined mass- 
and-energy spectrum: 

• Hadronic-mass moments are used as presented in 
this paper. We select the following subset for the 
fit: (m 2 x ) for p\ > 0.9, 1.1, 1.3, 1.5 GeV/c, (m 4 x ) 
forp* > 0.8,1.0,1.2,1.4GeV/c, and (m%) for p\ > 
0.9,l.l,1.3,1.5GeV/c. 

• Moments of the combined mass-and-energy spec- 
trum as presented in this paper. The follow- 
ing subset of moments is included in the fit: 
(n 2 x ) for p\ > 0.9,l.l,1.3,1.5GeV/c, (n\) for 
p\ > 0.8, 1.0, 1.2,1.4GeV/c, and (n x ) for p} > 
0.9,l.l,1.3,1.5GeV/c. 

Both fits include the updated lepton-energy moments as 
presented in this paper with radiative corrections as well 
as photon-energy moments measured in B — > X s j decays 
as presented in [14-16]. We use the partial branching 
fraction B p * measured for p* e > 0.6, 1.0, 1.5GcV/c and 

the moments'"^) for V* t > 0.6, 0.8, 1.0, 1.2, 1.5 GcV/c. 
The lepton-energy moments (Ef) are used for the min- 
imum lepton momentum p\ > 0.6, 1.0, 1.5 GeV/c and 
(Ef) for p\ > 0.8, 1.2GeV/c. We include the photon- 
energy moments (E^) for the minimum photon energies 
E 1 > 1.9 GeV and E 1 > 2.0 GeV, and (E*) for E 1 > 
1.9 GeV. In addition, we use tb = /o^o + (1 — fo)r± = 
(1.585 ± 0.007) ps, taking into account the lifetimes [21] 
of neutral and charged B mesons, To and t±, and their 
relative production rates, f = 0.491 ± 0.007 [21]. 

C. Theoretical Uncertainties 

As discussed in [17] and specified in [20], the following 
theoretical uncertainties are taken into account: 

The uncertainty related to the uncalculatcd perturba- 
tive corrections to the Wilson coefficients of nonperturba- 
tive operators are estimated by varying the corresponding 
parameters p 2 and by 20% and and p 3 LS by 30% 
around their expected values. Uncertainties for the per- 
turbative corrections are estimated by varying a s up and 
down by 0.1 for the hadronic-mass moments and by 0.04 
for the lepton-energy moments around its nominal value 
of a s = 0.22. Uncertainties in the perturbative correc- 
tions of the quark masses m b and m c are addressed by 
varying both by 20 MeV/c 2 up and down around their 
expected values. 

For the extracted value of | V^j,| an additional error of 
1.4% is added for the uncertainty in the expansion of the 



scmilcptonic rate Tsl [11, 55]. It accounts for remain- 
ing uncertainties in the perturbative corrections to the 
leading operator, uncalculated perturbative corrections 
to the chromomagnetic and Darwin operator, higher or- 
der power corrections, and possible nonperturbative ef- 
fects in the operators with charm fields. This uncertainty 
is not included in the theoretical covariance matrix Chqe 
but is listed separately as a theoretical uncertainty on 
\V cb \. 

For the predicted photon-energy moments (E™), addi- 
tional uncertainties are taken into account. As outlined 
in [52], uncertainties of 30% of the applied bias correction 
to the photon-energy moments and half the difference 
in the moments derived from two different distribution- 
function ansatze have to be considered. Both contribu- 
tions are added linearly [17]. 

The theoretical covariance matrix Chqe is constructed 
by assuming fully correlated theoretical uncertainties 
for a given moment with different lepton-momentum or 
photon-energy cutoffs and assuming uncorrelated theo- 
retical uncertainties for moments of different orders and 
types. The additional uncertainties considered for the 
photon-energy moments arc assumed to be uncorrelated 
for different moments and photon-energy cutoffs. 

D. Results 

In the following, the results of the two fits, one includ- 
ing the measurement of hadronic-mass moments and the 
other including the measured moments of the combined 
mass-and-energy spectrum instead, are discussed. 

We use a parameterized MC simulation to separate fit 
parameter uncertainties into experimental and theoreti- 
cal contributions. The simulation uses a set of expected 
moments randomly varied with either Ctot or C eX p- Fits 
to these moments allow for the determination of the ex- 
pected total and experimental uncertainties, respectively. 
The final experimental and theoretical uncertainties are 
calculated from the final total uncertainties by means of 
their simulated relative expected fractions. 



1. Combined Fit Including Hadronic-Mass Moments 

A comparison of the fit including hadronic-mass mo- 
ments with the measured moments is shown in Fig. 12. 
The moments (m x ) and (m x ) as well as the com- 
bined mass-and-energy moments are not included in the 
fit and thus provide an unbiased comparison with the fit- 
ted HQE prediction. We find an overall good agreement, 
also indicated by \ 2 = 10.9 for 28 degrees of freedom. 
Results for the SM and HQE parameters extracted from 
the fit are summarized in Table III. We find |U c fc| = 
(42.05 ±0.45 ±0.70) x 10" 3 , B{B -> X c e~D) = (10.64 ± 
0.17 ±0.06)%, m b = (4.549 ± 0.031 ± 0.038) GeV/c 2 , 
and m c = (1.077 ± 0.041 ± 0.062) GcV/c 2 , where the er- 
rors correspond to experimental and theoretical uncer- 
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taintics, respectively. The fitted quark masses have a 
large correlation of 95% resulting in a more precise de- 
termination of the quark mass difference, m b — m c = 
(3.472 ± 0.032) GeV/c 2 , where the error is the total un- 
certainty. We translate the quark masses which were ex- 
tracted in the kinetic scheme into the MS scheme us- 
ing calculations up to 0{a 2 s ) accuracy [11]. The transla- 
tion yields m b (m b ) = (4. 186 ± 0.044 ± 0.015) GcV/c 2 and 
m c (m c ) = (1.196 ±0.059 ±0.050) GcV/c 2 , where the first 
uncertainty is a translation of the uncertainty obtained 
in the kinetic scheme and the second corresponds to an 
estimate for the uncertainty of the transformation itself. 



2. Combined Fit Including Combined Mass-and-Energy 
Moments 

Figure 13 shows a comparison of the measured mo- 
ments and the fit including the measured combined mass- 
and-energy moments. We find an overall good agree- 
ment with x 2 = 8.2 for 28 degrees of freedom. The 
fit yields predictions of the hadronic-mass moments that 
are in good agreement with the measurement. Nu- 
merical results of the fit are summarized in Table IV. 
We find \Vcb\ = (41.91 ± 0.48 ± 0.70) x 10~ 3 , B(B -> 
X c e~v) = (10.64±0.17±0.06)%, m b = (4.566 ± 0.034 ± 
0.041) GeV/c 2 , and m c = (1.101 ± 0.045 ± 0.064) GcV/c 2 , 
where the errors correspond to experimental and theo- 
retical uncertainties, respectively. The two masses are 
correlated with 95%. Their difference is m b — m c = 
(3.465 ± 0.032) GeV/c 2 , where the stated uncertainty 
corresponds to the total uncertainty. The extracted 
quark masses translate into the MS scheme using [11] as 
m b (m b ) = (4.201 ± 0.047 ± 0.015) GcV/c 2 and m c (m c ) = 
(1.215±0.062±0.050) GcV/c 2 , where the first uncertainty 
is a translation of the uncertainty obtained in the kinetic 
scheme and the second corresponds to an estimate for 
the uncertainty of the transformation itself. 



3. Comparison of Results 

Comparing the result of the fit that includes moments 
of the n 2 x distribution with that including hadronic-mass 
moments instead, we find good agreement of all fit pa- 
rameters and their uncertainties. The differences be- 
tween the fit values are 0.2 a for \V c b\, 0.3 a for m b , 
and 0.3 a for m c . The uncertainties of all fit parameters 
in both fits agree within 8%. 

Figure 14 shows A\ 2 — 1 contours of both fits in the 
(m b , \ V cb \) and (to;,,/! 2 ) planes. We find an almost iden- 
tical precision for the fit values of \V cb \, m b , and /i 2 . In 
the Figure, we also show the results of two fits with re- 
duced sets of input measurements. To illustrate the influ- 
ence of the photon-energy measurements, a fit with only 
hadronic-mass and lepton-encrgy moments is performed. 
For further comparison we also perform a fit with only 
hadronic-mass moments and partial branching fractions. 



The fits with reduced experimental input show a signifi- 
cantly reduced accuracy of the extracted parameters. 

As our primary results we choose the values extracted 
from the fit with hadronic-mass moments since this fit 
has been used extensively before. Its results are in good 
agreement with earlier determinations [17, 59], but their 
uncertainties arc slightly larger because of the restrictions 
to BABAR data only. 

The use of combined mass-and-energy moments n 2 x 
does not lead to a more precise determination of the fun- 
damental physics parameters |V c fc|, m b , and m c . How- 
ever, the agreement of both fits confirms that higher- 
order corrections, which are needed for the expansion of 
the hadronic-mass moments but not for the n 2 x moments, 
have been estimated correctly. A significant change in the 
uncertainties of the SM and HQE parameters would have 
indicated a too naive treatment of the corrections for the 
mass moments [57]. Consequently, the presented results 
have increased the confidence into the validity of error es- 
timates that have to be made for a reliable determination 
of m b , m c , and |V^&|- 
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FIG. 12: The measured hadronic-mass moments (m x ), combined mass-and-energy moments (rix), electron-energy moments 
(Ei), partial branching fractions B, and photon-energy moments {E"), as a function of the minimum lepton momentum p| min 
and minimum photon energy E Jtrnin compared with the result of the simultaneous fit (solid line) to moments of the hadronic 
mass spectrum, electron-energy moments, and photon-energy moments. The solid data points mark the measurements included 
in the fit. Moments of semileptonic decays B — > X c l~V are marked by (•). Photon-energy moments of Ref. [14] are marked 
by (■), of Ref. [15] by (•), and of Ref. [16] by (*). Open data points are not used in the fit. The vertical bars indicate the 
experimental errors. The dashed lines correspond to the total fit uncertainty as obtained by converting the fit errors of each 
individual HQE parameter into an error for the individual moment. 
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FIG. 13: The measured hadronic-mass moments (my, combined mass-and-energy moments (n^), electron-energy moments 
(-B*), partial branching fractions B, and photon-energy moments (E™), as a function of the minimum lepton momentum 
Pl.min an d minimum photon energy E^, m in compared with the result of the simultaneous fit (solid line) to moments of the 
combined mass-and-energy spectrum, electron-energy moments, and photon-energy moments. The solid data points mark the 
measurements included in the fit. Moments of semileptonic decays B — * X C £~V are marked by (•). Photon-energy moments of 
Ref. [14] are marked by (■), of Ref. [15] by (•), and of Ref. [16] by (~k)- Open data points are not used in the fit. The vertical 
bars indicate the experimental errors. The dashed lines correspond to the total fit uncertainty as obtained by converting the 
fit errors of each individual HQE parameter into an error for the individual moment. 



23 



TABLE III: Results of the simultaneous fit to moments of the hadronic-mass spectrum, electron-energy moments, and photon- 
energy moments. For \V cb \ we account for an additional theoretical uncertainty of 1.4% from the uncertainty in the expansion 
of Tsl- Correlation coefficients for all parameters are summarized below the central values. 
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TABLE IV: Results of the simultaneous fit to moments of the combined mass-and-energy spectrum, electron-energy moments, 
and photon-energy moments. For |V c t,| we account for an additional theoretical uncertainty of 1.4% from the uncertainty in the 
expansion of Tsl- Correlation coefficients for all parameters are summarized below the central values. 





x 10 3 


m b [GcV/c"] 


m c [GeV/c"] 


b[%] pi 


[GeV*] fx' G 


[GeV*] p A D 


[GeV d ] pi s 


[GeV J ] 


Results 


41.91 


4.566 


1.101 


10.637 


0.452 


0.304 


0.190 


-0.156 


A ex p 


0.48 


0.034 


0.045 


0.166 


0.023 


0.047 


0.013 


0.079 


Atheo 


0.38 


0.041 


0.064 


0.061 


0.065 


0.039 


0.031 


0.052 


Ar sl 


0.59 
















Atot 


0.85 


0.053 


0.078 


0.176 


0.069 


0.061 


0.034 


0.095 


\v cb \ 


1.00 


-0.43 


-0.24 


0.74 


0.41 


-0.43 


0.43 


0.15 


m b 




1.00 


0.95 


0.04 


-0.58 


0.20 


-0.30 


-0.27 


m c 






1.00 


0.11 


-0.62 


-0.05 


-0.30 


-0.19 


B 








1.00 


0.17 


-0.09 


0.16 


-0.05 


2 

/'it 










1.00 


0.01 


0.68 


0.14 














1.00 


-0.05 


-0.05 


pi 














1.00 


-0.08 


pis 
















1.00 



24 




FIG. 14: Ax 2 = 1 contours for different fits in the (m&, \V c b\) an d (mi,,/ij) planes. We compare the results of the two fits 
including the full sets of measured moments, one including hadronic-mass moments (black line) and one including moments of 
the rix distribution instead (blue dashed line) , with a fit including only hadronic-mass and lepton-energy moments (red dotted 
line) and a fit including only hadronic-mass moments and partial branching fraction measurements (magenta dashed-dotted 
line). We do not include the additional uncertainty of 1.4% due to the expansion of Tsl in the plotted values of \V c b\- 
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TABLE A.I: Results for the moments (mjj) with k = 1 . . . 3 for different minimum lepton momenta p* t min with statistical and 
systematic uncertainties. The systematic uncertainties are grouped in five categories having related sources: MC statistics 
contains the statistical uncertainties of the calibration curves and of the residual background. Simulation related is the sum 
of uncertainties due to neutral and charged reconstruction efficiency differences in data and MC, particle identification, and 
mismodeling of final state radiation. The category extraction method contains the conservative estimate of half of the bias 
correction. The category background sums all contributions from the variation of the residual background components. The 
category signal model sums the impact of the variation of the signal decay branching fractions. Minimum lepton momenta are 
given in GeV/c. Moments and uncertainties are given in (GeV/c 2 ) fe . 
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1.5 


2.0494 


±0.0081 


±0.0198 


0.0036 


0.0074 


0.0168 


0.0061 


0.0019 




1.6 


2.0430 


±0.0092 


±0.0221 


0.0038 


0.0082 


0.0187 


0.0070 


0.0018 




1.7 


2.0387 


±0.0109 


±0.0265 


0.0047 


0.0081 


0.0232 


0.0083 


0.0015 




1.8 


2.0370 


±0.0143 


±0.0337 


0.0069 


0.0097 


0.0299 


0.0098 


0.0013 




1.9 


2.0388 


±0.0198 


±0.0413 


0.0082 


0.0123 


0.0355 


0.0150 


0.0008 


2 


0.8 


4.429 


±0.029 


±0.070 


0.027 


0.047 


0.030 


0.018 


0.008 




0.9 


4.416 


±0.027 


±0.063 


0.020 


0.041 


0.033 


0.016 


0.008 




1.0 


4.394 


±0.026 


±0.058 


0.020 


0.033 


0.035 


0.015 


0.008 




1.1 


4.354 


±0.026 


±0.063 


0.019 


0.031 


0.043 


0.016 


0.008 




1.2 


4.308 


±0.026 


±0.058 


0.019 


0.030 


0.039 


0.015 


0.007 




1.3 


4.281 


±0.027 


±0.061 


0.020 


0.029 


0.044 


0.016 


0.007 




1.4 


4.253 


±0.028 


±0.066 


0.021 


0.028 


0.051 


0.018 


0.006 




1.5 


4.220 


±0.031 


±0.070 


0.015 


0.029 


0.058 


0.019 


0.006 




1.6 


4.183 


±0.034 


±0.078 


0.015 


0.032 


0.065 


0.022 


0.005 




1.7 


4.158 


±0.040 


±0.094 


0.019 


0.032 


0.082 


0.026 


0.004 




1.8 


4.145 


±0.051 


±0.120 


0.026 


0.036 


0.107 


0.031 


0.004 




1.9 


4.136 


±0.069 


±0.142 


0.031 


0.046 


0.122 


0.048 


0.002 


3 


0.8 


9.57 


±0.11 


±0.25 


0.11 


0.18 


0.07 


0.06 


0.02 




0.9 


9.49 


±0.10 


±0.22 


0.08 


0.15 


0.08 


0.05 


0.03 




1.0 


9.41 


±0.09 


±0.18 


0.06 


0.11 


0.09 


0.04 


0.03 




1.1 


9.25 


±0.09 


±0.19 


0.06 


0.10 


0.11 


0.04 


0.03 




1.2 


9.09 


±0.08 


±0.18 


0.06 


0.10 


0.10 


0.04 


0.02 




1.3 


8.98 


±0.08 


±0.18 


0.06 


0.09 


0.11 


0.04 


0.02 




1.4 


8.88 


±0.09 


±0.19 


0.06 


0.09 


0.13 


0.04 


0.02 




1.5 


8.75 


±0.09 


±0.19 


0.04 


0.09 


0.15 


0.05 


0.02 




1.6 


8.61 


±0.10 


±0.22 


0.05 


0.10 


0.17 


0.06 


0.02 




1.7 


8.51 


±0.11 


±0.26 


0.05 


0.10 


0.22 


0.07 


0.01 




1.8 


8.45 


±0.14 


±0.33 


0.07 


0.11 


0.29 


0.08 


0.01 




1.9 


8.38 


±0.19 


±0.38 


0.08 


0.13 


0.32 


0.12 


0.00 
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TABLE A. II: Results for the moments (my with k = 4. . . 6 for different minimum lepton momenta p*i imin with statistical and 
systematic uncertainties. The systematic uncertainties are grouped in five categories having related sources: MC statistics 
contains the statistical uncertainties of the calibration curves and of the residual background. Simulation related is the sum 
of uncertainties due to neutral and charged reconstruction efficiency differences in data and MC, particle identification, and 
mismodeling of final state radiation. The category extraction method contains the conservative estimate of half of the bias 
correction. The category background sums all contributions from the variation of the residual background components. The 
category signal model sums the impact of the variation of the signal decay branching fractions, moment measurements. 
Minimum lepton momenta are given in GeV/c. Moments and uncertainties are given in (GeV/c 2 ) fe . 



k 




(m«x) 




&sys 


MC 


simulation 


extraction 


back- 


signal 




[GeV/c ] 








statistics 


related 


method 


groud 


model 


4 


0.8 


21.20 


±0.39 


±0.84 


0.35 


0.61 


0.14 


0.19 


0.11 




0.9 


20.83 


±0.33 


±0.69 


0.26 


0.51 


0.17 


0.15 


0.11 




1.0 


20.55 


±0.30 


±0.56 


0.24 


0.35 


0.19 


0.12 


0.12 




1.1 


20.01 


±0.27 


±0.55 


0.19 


0.32 


0.27 


0.11 


0.12 




1.2 


19.48 


±0.25 


±0.49 


0.17 


0.29 


0.23 


0.09 


0.10 




1.3 


19.09 


±0.25 


±0.52 


0.17 


0.33 


0.27 


0.10 


0.07 




1.4 


18.77 


±0.25 


±0.52 


0.17 


0.29 


0.32 


0.11 


0.07 




1.5 


18.33 


±0.26 


±0.50 


0.13 


0.24 


0.37 


0.11 


0.06 




1.6 


17.85 


±0.27 


±0.55 


0.12 


0.27 


0.42 


0.13 


0.05 




1.7 


17.50 


±0.30 


±0.66 


0.14 


0.26 


0.56 


0.15 


0.03 




1.8 


17.28 


±0.37 


±0.83 


0.18 


0.27 


0.73 


0.18 


0.03 




1.9 


16.99 


±0.48 


±0.90 


0.21 


0.34 


0.76 


0.27 


0.01 


5 


0.8 


48.51 


±1.39 


±2.90 


1.37 


2.10 


0.15 


0.64 


0.51 




0.9 


46.87 


±1.14 


±2.21 


0.84 


1.67 


0.31 


0.46 


0.49 




1.0 


46.00 


±0.97 


±1.74 


0.79 


1.07 


0.36 


0.32 


0.50 




1.1 


44.20 


±0.85 


±1.61 


0.57 


0.94 


0.60 


0.30 


0.48 




1.2 


42.55 


±0.77 


±1.44 


0.53 


0.88 


0.50 


0.24 


0.37 




1.3 


41.29 


±0.72 


±1.47 


0.43 


1.01 


0.61 


0.24 


0.28 




1.4 


40.31 


±0.70 


±1.45 


0.47 


0.94 


0.74 


0.25 


0.26 




1.5 


38.88 


±0.70 


±1.26 


0.33 


0.65 


0.84 


0.26 


0.23 




1.6 


37.35 


±0.70 


±1.38 


0.34 


0.71 


1.00 


0.29 


0.15 




1.7 


36.28 


±0.78 


±1.61 


0.34 


0.68 


1.32 


0.34 


0.10 




1.8 


35.56 


±0.94 


±2.00 


0.47 


0.69 


1.73 


0.41 


0.08 




1.9 


34.58 


±1.18 


±2.11 


0.56 


0.86 


1.73 


0.61 


0.04 


6 


0.8 


115.20 


±4.73 


±11.43 


4.39 


6.84 


5.64 


2.02 


3.76 




0.9 


107.97 


±3.74 


±8.32 


2.54 


5.36 


3.74 


1.36 


3.22 




1.0 


105.19 


±3.09 


±6.19 


2.34 


3.27 


2.26 


0.87 


3.05 




1.1 


99.35 


±2.60 


±5.19 


1.90 


2.85 


0.81 


0.79 


2.80 




1.2 


94.82 


±2.28 


±4.35 


1.53 


2.49 


0.23 


0.57 


2.16 




1.3 


91.01 


±2.05 


±4.09 


1.36 


2.73 


0.31 


0.53 


1.73 




1.4 


88.02 


±1.94 


±3.86 


1.23 


2.57 


0.94 


0.53 


1.57 




1.5 


83.46 


±1.86 


±3.35 


0.88 


2.20 


1.07 


0.53 


1.47 




1.6 


78.84 


±1.81 


±3.17 


0.84 


1.85 


1.73 


0.61 


0.97 




1.7 


75.87 


±1.98 


±3.92 


0.91 


1.73 


3.10 


0.76 


0.29 




1.8 


73.66 


±2.35 


±4.70 


1.06 


1.69 


4.05 


0.91 


0.25 




1.9 


70.70 


±2.83 


±4.77 


1.23 


2.09 


3.86 


1.33 


0.14 
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TABLE A. Ill: Results for (njj) for k = 2,4,6 for all minimum lepton momentum values p| m i„. The statistical uncertainty 
contains the uncertainty arising from the limited data sample and an additional statistical uncertainty arising from the 
determination of the combinatorial background. The systematic uncertainties are grouped in five categories having related 
sources: MC statistics contains the statistical uncertainties of the calibration curves and of the residual background. Simulation 
related is the sum of neutral and charged reconstruction efficiency differences in data and MC, J3 m i ss — cp m i ss differences, 
mismodeling of final state radiation, and PID impact. The category extraction method contains the conservative estimate of 
half of the bias correction and the impact of the calibration curve binning. The category background sums all contributions 
from the variation of the residual background components. The category signal model sums the impact of the variation of the 
signal decay branching fractions. 



k Pl.min 


(n k x) 




0"stat 




<^sys 


MC 


simulation extraction 


back- 


signal 


[GeV/c] 












statistics 


related 


method 


groud model 


2 0.8 


1.522 


± 


0.049 


± 


0.056 


0.020 


0.050 


0.011 


0.012 


0.004 


0.9 


1.483 


± 


0.047 


± 


0.057 


0.015 


0.054 


0.009 


0.009 


0.004 


1.0 


1.465 


± 


0.044 


± 


0.041 


0.013 


0.037 


0.009 


0.008 


0.003 


1.1 


1.438 


± 


0.037 


± 


0.040 


0.012 


0.037 


0.009 


0.006 


0.003 


1.2 


1.449 


± 


0.034 


± 


0.038 


0.011 


0.036 


0.006 


0.005 


0.002 


1.3 


1.428 


± 


0.031 


± 


0.031 


0.010 


0.027 


0.006 


0.006 


0.004 


1.4 


1.400 


± 


0.030 


± 


0.028 


0.009 


0.025 


0.006 


0.006 


0.004 


1.5 


1.369 


± 


0.035 


± 


0.032 


0.009 


0.029 


0.008 


0.007 


0.005 


1.6 


1.346 


± 


0.033 


± 


0.027 


0.009 


0.020 


0.014 


0.007 


0.003 


1.7 


1.344 


± 


0.037 ± 


0.029 


0.010 


0.020 


0.015 


0.008 


0.004 


1.8 


1.337 


± 


0.038 


± 


0.035 


0.013 


0.014 


0.029 


0.008 


0.004 


1.9 


1.196 


± 


0.032 


± 


0.033 


0.017 


0.017 


0.020 


0.009 


0.003 


4 0.8 


3.54 


± 


0.41 


± 


0.39 


0.14 


0.34 


0.08 


0.10 


0.03 


0.9 


3.21 


± 


0.37 


± 


0.36 


0.11 


0.32 


0.09 


0.05 


0.02 


1.0 


3.00 


± 


0.29 


± 


0.25 


0.09 


0.21 


0.09 


0.04 


0.02 


1.1 


2.74 


± 


0.22 


± 


0.17 


0.06 


0.14 


0.09 


0.01 


0.02 


1.2 


2.81 


± 


0.19 


± 


0.20 


0.06 


0.15 


0.12 


0.02 


0.03 


1.3 


2.60 


± 


0.15 


± 


0.16 


0.04 


0.10 


0.11 


0.01 


0.04 


1.4 


2.51 


± 


0.13 


± 


0.12 


0.04 


0.08 


0.09 


0.01 


0.03 


1.5 


2.34 


± 


0.13 


± 


0.13 


0.03 


0.09 


0.09 


0.00 


0.02 


1.6 


2.11 


± 


0.10 


± 


0.09 


0.03 


0.06 


0.06 


0.00 


0.01 


1.7 


2.03 


± 


0.12 


± 


0.08 


0.03 


0.06 


0.04 


0.00 


0.01 


1.8 


1.98 


± 


0.10 


± 


0.06 


0.04 


0.04 


0.02 


0.00 


0.01 


1.9 


1.57 


± 


0.07 


± 


0.06 


0.03 


0.05 


0.02 


0.00 


0.01 


6 0.8 


13.52 


± 


3.93 


± 


3.42 


1.37 


2.97 


0.49 


0.81 


0.34 


0.9 


10.87 


± 


2.78 


± 


2.65 


0.93 


2.39 


0.52 


0.37 


0.24 


1.0 


9.02 


± 


2.22 


± 


1.88 


0.84 


1.55 


0.54 


0.29 


0.20 


1.1 


7.06 


± 


1.35 


± 


0.78 


0.35 


0.58 


0.34 


0.07 


0.14 


1.2 


7.50 


± 


1.16 


± 


0.92 


0.32 


0.68 


0.49 


0.11 


0.18 


1.3 


6.28 


± 


0.84 


± 


0.64 


0.22 


0.38 


0.41 


0.06 


0.20 


1.4 


5.83 


± 


0.62 


± 


0.49 


0.16 


0.30 


0.32 


0.06 


0.12 


1.5 


4.99 


± 


0.49 


± 


0.52 


0.13 


0.30 


0.40 


0.03 


0.05 


1.6 


3.93 


± 


0.32 


± 


0.31 


0.11 


0.16 


0.24 


0.03 


0.03 


1.7 


3.63 


± 


0.35 


± 


0.22 


0.09 


0.13 


0.15 


0.02 


0.03 


1.8 


3.42 


± 


0.23 


± 


0.19 


0.09 


0.11 


0.12 


0.02 


0.03 


1.9 


2.51 


± 


0.16 


± 


0.13 


0.06 


0.11 


0.03 


0.02 


0.02 



